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ABSTRACT 

One  of  the  non-linearities  which  occur  in  feedback  control  systems 
is  transportation  lag.  All  of  the  currently  popular  mathematical  tech- 
niques available  for  the  analysis  of  these  systems  are  severely  compli- 
cated by  the  introduction  of  transportation  lag.   The  suitability  of  a 
numerical  method  employing  the  Z- trans  format  ion  to  this  problem  was  in- 
vestigated, as  well  as  the  particular  advantages  and  disadvantages  of 
the  other  methods.   Solutions  by  the  numerical  method  were  compared 
with  solutions  by  other  methods  and  the  results  obtained  presented  in 
tabular  form.   These  results  demonstrate  the  superiority  of  the  numeri- 
cal method  in  determining  the  transient  response  of  these  systems  with- 
in reasonable  engineering  accuracy  in  a  much  shorter  time. 

The  writers  wish  to  express  their  appreciation  to  Dr.  M.  P.  Pastel 
for  his  guidance  during  the  initial  stages  of  the  work,  and  to  Dr.  G.  J. 
Thaler  for  his  advice  and  suggestions  throughout  the  investigation  and 
his  aid  in  the  preparation  of  this  paper. 
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1.    Introduction 

1. 1   Purpose  of  the  Study 

Feedback  control  systems  are  seldom  truly  linear,  but  rather 
they  are  subject  to  the  introduction  of  a  considerable  number  of  non- 
linearities.   Some  non-linearities  may  be  introduced  intentionally  by 
the  system  designer.   An  example  of  this  might  be  the  selection  of  a 
relay  servomechanism  with  its  attendant  non-linearity,  when  a  linear 
system  could  accomplish  the  task  just  as  well.   On  the  other  hand  there 
may  be  no  way  of  avoiding  a  non-linearity  in  the  system  design.   An  ex- 
ample of  this  would  be  the  mandatory  inclusion  of  a  gear  train  with  its 
attendant  backlash. 

A  non-linearity  which  frequently  can  neither  be  avoided  nor 
neglected  is  transportation  lag. 

Because  of  the  variable  influence  of  transportation  lag  in 
feedback  control  systems  and  its  avoidance  whenever  possible,  it  has 
received  relatively  little  prominence  in  the  literature.   Rather,  the 
mathematical  tools  capable  of  analyzing  systems  with  transportation  lag 
have  been  indicated  without  further  amplification.   While  it  is  true  that 
the  tools  indicated  for  the  analysis  of  systems  with  transportation  lag 
will  yield  satisfactory  results,  the  problem  does  not  end  there.   These 
standard  tools,  while  yielding  accurate  results  with  a  reasonable  amount 
of  work  for  linear  systems,  require  considerably  more  work  and  time  when 
the  problem  is  complicated  by  transportation  lag.   Consequently,  a  need 
for  a  better  tool  than  those  currently  in  use  for  the  solution  of  pro- 
blems involving  transportation  lag  is  indicated. 

In  1954,  a  paper  by  Ragazzini  and  Bergen  outlined  a  numerical 


methods  technique  for  determining  the  transient  response  of  linear  feed- 
back control  systems.   The  nature  of  the  mathematics  which  the  method 
used  indicated  that  it  might  handle  the  non-linearity  of  transportation 
lag  more  conveniently  than  the  aforementioned  standard  mathematical 
methods  generally  known  to  engineers.   The  purpose  of  this  study  is  to 
evaluate  the  numerical  methods  technique  of  Ragazzini  and  Bergen  in 
analyzing  feedback  control  systems  with  transportation  lag  and  to  com- 
pare the  accuracy  attainable  and  the  work  involved  in  using  this  method 
with  the  accuracy  attainable  and  work  involved  in  using  more  popular 
methods. 

In  doing  so,  it  will  be  necessary  to  describe  briefly  the 
common  analytical  and  graphical  methods  available  for  solving  linear  feed- 
back systems.   The  increased  difficulty  and  complexity  introduced  by  the 
transportation  lag  into  these  methods  will  be  demonstrated.   Finally, 
the  application  of  the  numerical  method  to  the  same  problem  will  be 
presented. 
1.2   Summary  of  the  Results  and  Conclusions 

The  numerical  methods  technique  was  compared  with  the  root 
locus,  Bode  diagram,  power  series  expansion,  and  second  order  Pade  approxi- 
mation approaches  to  the  problem  of  system  analysis.   The  Bode  diagram 
proved  to  be  the  speediest  method  of  all  for  determining  system  stability, 
but  this  method  yields  no  other  specific  information  regarding  the  over- 
all system  performance.   For  the  full  transient  response  the  numerical 
methods  technique  proved  to  be  faster  and  more  accurate  than  the  second 
order  Pade  approximation.   The  root  locus  method,  while  retaining  its 
potential  for  good  accuracy,  required  a  prohibitive  amount  of  time  and 
tedium  to  be  considered  a  practical  working  tool  for  problems  involving 
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transportation  lag.   The  power  series  expansion  approach  to  the  pro- 
blem proved  to  be  totally  impractical  for  other  than  extremely  simple 
systems. 

In  short,  the  numerical  methods  technique  is  the  fastest  and 
most  accurate  means  of  determining  the  transient  response  of  a  feedback 
control  system  with  transportation  lag.   Moreover,  the  method  can  be 
understood  and  applied  by  engineers  and  students  whose  background  includes 
only  a  first  course  in  linear  servomechanisms.   In  the  solution  of  the 
problem  by  the  numerical  methods  technique,  no  more  complicated  computing 
device  than  the  desk  calculator  is  required. 

In  all  cases,  solutions  determined  by  the  various  methods  were 
compared  with  the  transient  response  obtained  from  programming  the  Fourth 
Order  Pade  Approximation  of  the  Laplace  shift  operator  on  a  ten  channel 
Donner  analog  computer.   This  result  was  the  most  correct  obtainable 
under  the  conditions  of  time  and  equipment  available  for  the  study. 


2.    Description  of  the  System 

The  feedback  control  system  to  be  studied  is  shown  in  block 
form  in  figure  2-1.   In  order  to  carry  out  the  goals  of  this  thesis,  a 
number  of  requirements  were  set  for  this  system. 

A  frequency  response  diagram  crossover  of  one  radian  per 
second  for  all  variations  of  the  system  was  desired.   For  this  reason, 
the  value  of  K  in  block  one  will  always  be  equal  to  one. 

The  forward  loop  transfer  function  must  always  be  low  pass  in 
nature.   For  this  reason  the  following  values  of  G(s)  were  selected  in 
block  two.   These  values  also  give  us  convenient  first,  second,  and 
third  order  systems  to  investigate. 

G(s)  =  1  2.1 

G(s)  =   - 2.2 

0-/S  +  I 

G(s)  =   I 2.3 

(o.\s+i)(o.oSs+\) 

Block  three  contains  the  transportation  lag.   The  nature  of 
this  block  and  its  analytic  representation  will  be  discussed  in  detail  in 
section  4. 

The  system  will  use  unity  feedback. 
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3.    Linear  System  Analysis 

3. 1  Discussion 

Various  methods  of  analysis  and  solution  of  linear  feedback 
systems  are  available.   Two  graphical  methods  currently  in  use  supply 
the  necessary  information  required  for  the  design  and  analysis  of  these 
systems.   For  the  purposes  of  this  study,  the  transient,  or  time  varying 
response  characteristics,  are  more  important. 

3.2  Frequency  Response 

The  Bode  network  therorems  presented  a  new  approach  to  the 
standard  Nyquist  stability  criterion.   One  of  the  great  values  of  the 
Bode  diagram  is  the  simplicity  with  which  approximate  results  of  suf- 
ficient accuracy  may  be  obtained. 

The  Bode  diagram  for  the  linear  second  order  system  described 
in  section  2  is  shown  in  figure  3-1.   Phase  margin  and  gain  margin  values 
are  indicated,  and  reasonable  predictions  of  system  performance  may  be 
made  from  these  values.   This  graphical  presentation  of  the  open  loop 
frequency  response  is  plotted  knowing  the  location  of  the  open  loop  poles 
and  zeros. 

System  stability  may  be  inferred  from  the  values  of  gain  margin 
and  phase  margin  on  the  open  loop  frequency  response.   Information  from 
the  Bode  diagram  may  be  transferred  to  the  Nichols  chart,  and  the  closed 
loop  frequency  response  may  be  obtained.    Information  from  the  Nichols 
chart  should  agree  closely  with  that  from  a  basic  Nyquist  plot. 
3.3   Root  Locus 

A  graphical  method  exists  for  determining  the  locus  of  the 
roots  of  the  closed  loop  system  from  the  location  of  the  poles  and  zeros 
of  the  open  loop  system.   The  root  locus  method  essentially  demonstrates 
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the  transient  response  characteristics  of  the  system  rather  than  the 
frequency  response. 

The  root  locus  plot  for  the  linear  second  order  system  describ- 
ed in  section  2  is  shown  in  figure  3-2.   Since  this  method  indicates 
directly  the  effect  on  the  closed  loop  response  of  the  system  gain  from 
zero  to  infinity,  it  is  obvious  that  this  system  is  stable  for  all  values 
of  K.   This  may  also  be  inferred  from  inspection  of  figure  3-1,  which 
shows  that  the  phase  angle  curve  always  remains  above  the  180  degree 
phase  angle  line. 
3.4  Transient  Response 

The  transient  response  of  a  control  system  contains  nearly  all 
the  parameters  required  by  the  engineer  for  the  evaluation  of  its  perfor- 
mance.  A  typical  transient  response  is  shown  in  figure  3-3  and  the  vital 
parameters  are  indicated.   From  a  study  of  these  values,  an  estimate  of 
the  worth  of  the  control  system  may  be  made,  and  the  need  for  compensa- 
tion may  be  determined. 

The  parameters  which  will  be  of  interest  are  the  following: 
Rise  Time  -  Length  of  time  for  the  system  to  first  reach  a  point  within 

57.  of  the  desired  value. 
Peak  Overshoot  -  Value  in  percent  by  which  the  system  overshoots  the  de- 
sired value. 
Time  to  Peak  Overshoot  -  Length  of  time  to  reach  peak  overshoot  point. 
Settling  Time  -  Length  of  time  for  the  system  to  settle  within  57.  of  the 
desired  value. 
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Figure  3-2.  Non-Dimensionalized  Root  Locus  of  Linear  Second  Order  System 

G(s)  = 

sCs+10; 


q.ndq.no 
10 


4.    Transportation  Lag 
4. 1  Discussion 

"Transportation  lag"  is  a  term  frequently  used  in  the  process 
industry  to  describe  a  pure  time  delay.   Other  terms  used  to  describe 
this  phenomenon  include  "dead  time",  "delay  time",  and  "true  time  lag". 
Figure  4-1  outlines  an  arrangement  which  might  be  used  to  control  the 
thickness  of  steel  plate  produced  by  a  rolling  mill.   A  sensor  indicates 
the  actual  plate  thickness  in  voltage  form.   The  difference  between  its 
output  and  another  voltage  representing  the  desired  thickness  forms  the 
error  voltage.   The  amplified  error  ultimately  positions  the  rolls.   The 
transfer  function  from  plate  thickness  at  the  sensor  X  through  the  ampli- 
fier, motor,  and  gear  train  to  roll  position,  hence  thickness  at  the  rolls 
Y,  differs  in  no  sense  from  ordinary  transfer  functions  considered  pre- 
viously.  From  Y  to  X  however  a  new  element  is  introduced  because  a  finite 
time  must  elapse  before  a  change  in  thickness  at  Y  reaches  the  point  of 
measurement  at  X,  d  units  of  length  away.   If  the  steel  plate  proceeds  to 
the  left  with  a  constant  velocity  v,  any  time  function  of  thickness  y(t) 
resulting  from  changes  in  roll  position  will  be  reproduced  at  the  sensor 
without  distortion  d/v  time  units  later.   It  is  precisely  this  time  quan- 
tum that  the  term  transportation  lag  describes. 

Transportation  lag  may  be  expressed  by  the  following  mathe- 
matical relationship:   Given  f(t)  for  t  ^0,  the  function  f(t-T)  differs 
from  the  function  f(t)  by  the  dead  time  T,  where  f(t-T)  is  defined  to  be 
zero  for  t^T. 

The  setup  of  problems  on  computers  occasionally  requires  the 
representation  of  transportation  lag  or  an  expression  of  the  form  f(t-T). 
The  brute- force  approach  to  the  solution  of  this  problem,  and  one  that 
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has  been  used  in  the  past,  is  to  generate  f(t)  in  the  computer  and  plot 
f(t)  vs.  t  on  a  specially  equipped  recording  device.   Schemes  involving 
such  mechanical  devices  to  accomplish  the  introduction  of  transportation 
lag  have  proved  to  be  workable  for  some  problems,  but  the  limitations 
imposed  by  the  physical  design  of  the  devices  and  the  introduction  of 
human  error  in  their  use  have  made  necessary  the  development  of  better 
methods. 

The  Laplace  transform  of  f(t-T)  is 

X[f(t-T)]  -  X[f(t)]c-sT=  F(S)CST 

-sT 
if  f(t-T)  =  0  for  t<T.   In  equation  4.1  e     is  commonly  called  the 

Laplace  shift  operator.   The  problem  of  generating  a  function  f(t-T) 

-sT 
can  therefore  be  reduced  to  the  problem  of  generating  the  function  e  '  . 

-sT 
Since  e    cannot  be  expressed  as  a  rational  fraction,  the  transfer  func- 
tion for  dead  time  can  be  simulated  only  approximately.   Thus  the  remain- 

-sT 
ing  methods  depend  on  an  approximation  of  the  transfer  function  e  '   by 

an  approximate  series  or  technique  of  network  synthesis. 

4.2  Methods  of  Simulation 

4.21  Power  Series 

The  best  known  approximation  to  the  function  e   is  the  Taylor 

series  expansion 

e*  =  /  +  x  +  JL  +A  +  •  •  +-fir  4.2 

2.!  3/  Y)[ 

This  expansion  is  convergent  for  all  values  of  x,  but  unfortunately 
the  rate  of  convergence  is  small  for  x  large.   For  this  reason  the 
Taylor  series  is  not  well  suited  to  the  generation  of  transportation 
lag  for  problems  involving  high  frequencies  or  long  time  constants. 
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-sT 
This  can  be  seen  by  replacing  s  by  ju)  in  the  expression  for  e  '  . 

Moreover,  in  attempting  to  determine  by  longhand  calculations,  the  trans- 
ient response  of  other  than  the  simplest  system  possible  involving  trans- 
portation lag,  the  introduction  of  the  Taylor  series  approximation  for 

-sT 
e  '   complicates  the  problem  so  severely  as  to  render  it  useless. 

4. 22  The  Pade  Approximation 

In  section  4.1  an  important  theorem  of  Laplace  transforms  was 

stated,  viz. : 

^[f  (t-T)j  =   X[-f(t)]e""S  repeated  4.1 

Thus  the  transfer  function  of  transportation  lag  is 

£Cffi-T)]        ST  Eo 

= ^    P     =    ■L-°  4  3 

l[Ut>]  L  Ei 

Substituting  s  =  jW  to  find  the  frequency  response  gives 


Thus 


]T    "    C*       =     cos  ooT  - frSiNLoT 


Ei 


4.4 


or  the  frequency  transfer  function  of  transportation  lag  is  characteriz- 
ed by  an  amplitude  ratio  which  is  unity  for  all  frequencies  and  a  phase 
angle  which  varies  linearly  with  frequency  over  the  desired  range  of  in- 
terest. 

-sT 
A  popular  method  of  approximating  the  transfer  function  e 

from  this  point  of  view  is  the  Pade  approximation.   It  stems  from  the 
fact  that  the  quotient  of  any  pair  of  complex-conjugate  polynomials  has 
an  amplitude  of  unity  and  also  possesses  considerable  phase  shift.   Ex- 
pressed in  its  basic  form,  the  Pade  approximation  is 
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where  ,   .  ^2. 

f     (y)   «   /  +   ^*    +    W"-/J  * 


4.6 


+ 


...+.   ^-')-  •  £•/,  x 4.7 

Gu,v(x)  =    /-      *4*f      +       afcv-Qx2- .... 

(v+u)/y    x     (v+Wf\/+u-<;  2/ 

.  ..  +  (-/]«   *(<*-»)  ■■-£-/  ^  4o8 

(V+u)  (V  +  u.-/^.  ••  (v+i)U! 

The  convergence  for  this  series  expansion  is  quite  rapid  and  often  values 

of  u  and  v  of  2  will  give  good  accuracy  for  short  time  lags  and  low  fre- 

-sT 
quencies.   The  Pade  approxitnant  of  e     for  u  =  v  =  2  is  given  in  equa- 
tion 4. 10  below. 


-sT  „     l-jSr(sT)    t  ^(ST') 
e    ~"     l+£(ST)    +^^T»)  4.9 

^>-sT^        sT2-fcs7--f-/£ 

v_     —   ~~"r 4.10 

In  terms  of  the  original  problem,  generating  f  (t-T)  when  f(t) 
is  known,  the  equation  for  the  second-order  Pade  approximant  can  be 
expressed  in  transfer- function  notation  as 

f(t-T)      _         D*T*-&DT-M£  4ai 

f(t)  D*T2+&DT+/2 

The  justification  of  the  validity  of  equation  4.11   is  demon- 
strated as  follows: 


L[j(t-T)]  =  F(i)c5T-Hm]c 


~$  repeated  4. 1 
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Substituting  equation  4.11   in  equation  4.1  gives 


4.12 


Rewriting  equation  4.12  gives 

(s2T'asT+/£)X[f(t-T)J  =  ($«?*- tsT+tt)-££«t>]         4.13 
Performing  the  inverse  Laplace  transformation  on  equation  4.13  gives 

(D*f+  (°DT+ 1 Z)f(t-T)  -  (D2T2-t>DT+IZ)f(t)  4.14 

Comparison  of  equations  4.14  and  4.12   reveals  that  the  equations  are 
identical.   This  is,  of  course,  a  logical  conclusion,  as  the  Laplacian 
operator  and  the  differential  operator  may  be  used  interchangeably  pro- 
viding the  initial  conditions  in  a  system  are  identically  zero.   This  is 

true  in  the  system  described  above. 

-sT 
The  fourth  order  Pade  approximant  (u  =  v  =  4)  for  e  '   is 

PsTc.         (ST?  -  ZOjsT)3  +   I  go  (sT)2-  g40  sT+  /4S0 
(sT)U2o(sTf+   I80(5T)V  MOsT+  /£>g0 

The  Pade  approximants  may  be  used  for  either  longhand  or  com- 
puter solutions.   For  use  in  longhand  solutions  the  approximant  to  be 

-sT 
used  is  merely  substituted  for  e   .   It  is  worthy  of  note  that  the 

order  of  the  characteristic  equation  is  raised  by  the  order  of  the 

-sT 
approximant  used  to  simulate  e   .   If  the  second  order  approximant  is 

used,  two  additional  roots  are  introduced  into  the  characteristic  equa- 
tion, and  if  the  fourth  order  approximant  is  used,  four  additional  roots 
are  introduced  into  the  characteristic  equation. 

The  technique  of  programming  the  Donner  Analog  Computer  for 
handling  the  transfer  function  of  the  fourth  order  Pade  approximant 
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is  included  in  Appendix  Section   9. 

The  method  of  simulating  transportation  lag  by  the  Pade  ap- 
proximants  yields  what  may  be  termed  optimally  flat  phase  characteris- 
tics at  the  origin  but  does  not  fit  the  phase  characteristics  as  well 
over  an  extended  frequency  range.   The  major  drawback  to  the  use  of  the 
Pade  approximants  to  simulate  transportation  lag  is  their  response  to  a 
step  input  or  other  input  having  high-frequency  components,  since  the 
amplitude  ratio  is  always  unity  but  the  phase  is  approximated  only  for 
lowcuT  products.   The  phase  shift  characteristics  of  certain  of  the  Pade 
approximants  are  shown  in  Figure  4-2.   As  a  result  of  the  inability  of 
these  approximants  to  linearly  phase  shift  high  frequencies,  the  output 
of  an  approximant  "rings"  badly  for  high-frequency  inputs.   This  behavior 
is  illustrated  in  Figure  4-3  which  depicts  the  response  of  the  fourth 
order  approximant  to  inputs  of  a  step  and  ramp  with  unity  feedback.   It 
is  noteworthy  that  the  "ringing"  is  considerably  less  for  the  ramp  in- 
put than  for  the  step  input.   The  reason,  of  course,  is  that  the  ramp  in- 
put contains  a  lower  percentage  of  high-frequency  terms  than  does  the  step 
input. 

In  the  solution  of  problems  concerning  feedback  control  systems 
with  transportation  lag,  advantage  may  be  taken  of  the  fact  that  the  typi- 
cal system  encountered  is,  in  effect,  a  low  pass  filter.   Block  diagram 
algebra  techniques  are  employed  in  the  solution  of  these  problems.   Since 
the  commutative  law  is  applicable  to  block  diagram  algebra,  the  actual 
location  of  the  transportation  lag  transfer  function  in  the  forward  path 
is  immaterial  and  it  may  be  placed  at  the  end  of  the  path  in  the  computer 
where  the  high  frequency  components  have  been  attenuated  to  the  greatest 
possible  extent.   Frequencies  which  have  been  attenuated  40  db.  or  more 
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CHART  NO.  RA- 2921-30  BRUSH      I IV  STRUM  E? 


Figure  4.-3 (a).  Response  of  Fourth  Order  Pade  Approximation  of 
0.9  Second  Time  Delay  to  Step  Input. 
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are  considered  to  be  negligible  and  whether  they  have  been  linearly 
phase  shifted  is  of  no  consequence.   This  tends  to  obviate  the  neces- 
sity of  having  to  linearly  phase  shift  a  high  vo   T  product.   Another 
consideration  which  makes  unnecessary  the  linear  phase  shift  of  high  uj  T 
products  in  the  solution  of  the  transient  response  of  feedback  control 
systems  is  that  high  uo  T  products  mean  extremely  large  phase  shifts  and 
probable  instability.   In  these  cases  the  transient  response  is  virtually 
meaningless. 

In  the  great  majority  of  cases  the  fourth  order  Pade  approxi- 
mant  is  more  than  adequate  for  the  simulation  of  transportation  lag.   In 
cases  where  linear  phase  shift  of  higher  ujT  products  than  those  shown  in 

figure  4-2  is  desirable,  higher  order  Pade  approximants  or  more  complex 

-sT 
curve  fitting  techniques  may  be  used  to  simulate  e  '  .   The  curve  fitting 

techniques  are  discussed  in  the  literature  on  analog  computers  but  the 
need  for  such  approaches  is  questionable.   These  techniques  were  devised 
because  it  was  felt  the  Pade  approximants  of  sixth  order  and  higher  fail- 
ed to  phase  shift  low  ojT  products  accurately.   In  the  course  of  this 
study  the  phase  shift  characteristics  of  the  sixth  and  eighth  order  Pade 
approximants  were  calculated.   The  results  are  shown  in  figure  4-4.   On 
the  basis  of  these  results  it  may  be  inferred  that  the  order  of  the  Pade 
approximation  may  be  increased  to  increase  the  magnitude  of  the  ^T  pro- 
duct that  is  linearly  phase  shifted  without  distorting  the  phase  shift  of 
low  tuT  products.   Because  the  Pade  approximation  may  be  expressed  in  a 
general  form  from  which  any  order  approximant  may  be  developed,  it  is 
felt  that  the  use  of  curve  fitting  techniques  requires  a  needless  expendi- 
ture of  effort  without  any  tangible  benefit. 
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5.    Effect  of  Time  Delay  on  Linear  System  Analysis 

The  addition  of  the  time  delay  into  the  linear  control  system 
under  investigation  introduces  several  factors  which  tend  to  make  the 
analysis  of  the  system  more  difficult.   Because  the  time  delay  may  basical- 
ly be  represented  by  an  infinite  series  of  terms,  this  creates  an  infinite 
number  of  roots.   The  effect  of  these  complications  on  the  standard  graph- 
ical   analysis  techniques  will  now  be  investigated.   In  all  cases,  the 
system  with  a  time  delay  of  1.0  seconds  will  be  demonstrated. 
5. 1  The  Frequency  Response 

The  standard  Bode  diagram  is  shown  in  figure  5-1  when  the  time 
delay  is  included  in  the  system.   The  difference  between  this  plot  and 
figure  3-1  is  readily  apparent.   The  phase  margin  has  been  reduced  by  a 
considerable  amount.   The  gain  margin  is  approximately  3.3  db  which 
compares  with  the  undelayed  system  which  was  stable  at  all  gains. 

The  time  delay  has  introduced  additional  branches  into  the  dia- 
gram.  These  branches  appear  due  to  the  periodic  nature  of  the  time  de- 
lay, and  the  increasing  phase  shift  introduced  at  higher  values  of  fre- 
quency.  It  is  evident  that  the  gain  margin  for  each  of  the  branches  in- 
creases.  This  is  analogous  to  the  statement  that  the  roots  of  the  pri- 
mary branch  would  go  unstable  first  if  the  gain  were  raised. 

The  construction  of  this  plot  is  simple  and  straightforward. 

For  the  system  in  question  the  open  loop  transfer  function  is 

-ST 

G0CS)  =   —^ 5.1 

s(o.is+i) 

The  gain  vs  logo;  curve  is  identical  to  the  undelayed  system 
curve  since  the  time  delay  introduces  no  additional  gain  into  the  system. 
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The  phase  angle  vs  logo)   curve  is  drastically  changed  since 
the  equation  for  the  phase  angle  is 

(f)   =   ^0°  +     *RC   rAfvJ  °'/oJ  +  w  "^  5.2 

It  is  the  presence  of  the  additional  term  of  the  frequency  times 
the  time  delay  which  introduces  the  additional  phase  shift  and  the  extra 
branches. 
5.2   Root  Locus 

The  identical  second  order  system  with  a  transportation  lag 
of  1.0  seconds  also  appears  greatly  changed  when  investigated  by  a 
root  locus  plot. 

Figure  3-2  showed  the  undelayed  second  order  system  which  is 
stable  for  all  gains.   Figure  5-2  has  the  additional  construction  de- 
tail which  we  will  require  to  obtain  the  root  locus.   In  the  same  manner 
as  the  minus  180  degree  phase  shift  line  was  obtained,  lines  of  constant 
phase  shift  of  170,  160,  150,  140  degrees,  and  so  forth,  may  be  plotted. 
In  addition,  lines  of  constant  phase  shift  caused  by  the  one  second  time 
delay  are  also  shown.   In  the  non-dimensionalized  plot  shown,  the  time  de- 
lay will  produce  a  complete  360  degree  phase  shift  every  0.628  radians. 

The  resulting  root  locus  is  shown  in  figure  5-3.   It  is  plot- 
ted by  combining  points  of  the  basic  system  and  the  time  delay  which 
total  to  a  180  degree  phase  shift.   It  is  evident  that  this  system  will 
go  unstable  with  only  a  small  increase  in  gain,  and  as  a  matter  of  fact, 
the  gain  at  the  stability  limit  just  before  the  locus  crosses  the  axis  is 
approximately  1.48,  the  same  value  of  gain  margin  obtained  from  the  Bode 
diagram. 
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Figure  5-2.  Root  Locus  Construction  Detail 
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An  extremely  interesting  result  of  this  plot  is  that  upper 
branches  of  the  root  locus  appear.   Only  two  of  these  are  shown  although 
an  infinite  number  exist.   Here  again  they  are  a  result  of  the  360  degree 
periodicity  introduced  by  the  time  lag.   It  is  also  noted  that  at  the 
same  gain,  the  additional  roots  on  the  upper  branches  are  attenuated  more 
than  the  root  on  the  primary  branch.   This  substantiates  the  Bode  diagram 
findings,  which  show  that  the  roots  on  the  primary  branch  would  go  unstable 
first. 

A  continuation  of  the  graphical  solution  described  above  for 
other  values  of  time  delay  will  yield  other  root  loci  curves  similar  to 
the  one  in  figure  5-3. 

The  gain  at  the  stability  limit  for  each  time  delay  was  calculat- 
ed.  The  product  of  the  time  constant  of  the  system  and  the  gain,  K  t^  , 
was  plotted  against  the  time  delay  on  logarithmic  cross-section  paper. 
The  result  is  a  straight  line  as  shown  in  figure  5-4.   It  was  therefore 
possible  to  write  a  single  equation  which  will  describe  the  gain  at  the 
stability  limit  in  terms  of  the  time  constant  of  the  second  order  system 
and  the  time  delay. 

La  lD 

Thus  for  second  order  systems  if  two  of  the  above  parameters  are 
known,  the  maximum  value  that  the  third  parameter  may  have  can  readily  be 
calculated. 
5.3  Transient  Response 

Neither  of  the  graphical  methods  discussed  here  will  directly 
produce  the  transient  response  of  any  system.   The  Bode  diagram  will  only 
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yield  values  from  which  an  approximate  transient  response  may  be  infer- 
red.  The  root  locus  method  provides  the  roots  of  the  system  on  the  s 
plane,  but  then  the  results  must  be  transformed  to  the  time  domain  and  the 
equation  of  the  output  determined.   The  transient  response  must  be  obtain- 
ed by  substituting  values  into  this  equation  and  obtaining  the  output  at 
each  specific  time.   Complicating  this  procedure  is  the  fact  that  an  in- 
finite number  of  roots  exist  on  the  root  locus,  and  enough  of  these  roots 
must  be  considered  to  obtain  reasonable  accuracy.   Each  time  an  addition- 
al set  of  these  complex  conjugate  roots  is  considered,  additional  time  and 
complexity  are  added  to  the  above  procedure. 

The  transient  response  of  this  delayed  second  order  system  is 
shown  in  figure  5-5.   It  is  here  that  the  influence  of  the  time  delay  on 
the  second  order  system  is  most  apparent.   The  original  second  order 
system  was  an  overdamped  system.   It  had  no  overshoot.   The  time  delay  has 
deteriorated  the  system  to  the  point  where  the  overshoot  is  approximately 
70  percent. 

It  is  interesting  to  note  that  a  time  delay  of  about  1.4  seconds 
would  cause  this  system  to  be  completely  unstable. 

The  effect  of  the  time  delay  on  the  transient  response  there- 
fore is  to  increase  overshoot  and  settling  time.   It  can  make  a  normally 
overdamped  system  highly  oscillatory  with  a  delay  of  only  one  second  in- 
troduced. 
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6.    Transient  Response  by  Numerical  Methods 

6. 1  Discussion 

A  mathematical  technique  now  exists  for  the  direct  solution  of 
the  output  of  a  linear  continuous  feedback  system.   An  extension  of  this 
numerical  method  to  the  solution  of  systems  with  transportation  lag  is 
easily  accomplished. 

6.2  The  Z  -  Transformation 

A  continuous  time  function  f(t)  is  sampled  by  means  of  a  sampl- 
ing switch  S  at  equal  sampling  intervals  of  time  T  in  figure  6-1.   The 
output  of  the  sampler  f*(t)  is  a  sequence  of  pulses  which  may  be  represent- 
ed for  mathematical  purposes  as  a  series  of  impulses  or  delta  functions 
whose  areas  are  equal  to  the  amplitude  of  f(t)  at  the  respective  sampl- 
ing instants.   Thus 

The  Laplace  transform  of  the  sequence  of  impulses  is  given  as 

I[f(t)hl4(«T)e-"T5 

10  =  0 

The  Laplace  transform  of  such  pulsed  functions,  therefore,  is  in 

~nTs 
terms  of  e    ,   Replacing  this  exponential  term  by  a  new  variable  Z,  where 

Z  =  e     accounts  for  the  name  Z- transform.   One  useful  characteristic 
of  the  Z-transform  is  that  the  infinite  series  in  equation  6.2  may  usual- 
ly be  expressed  in  closed  form  so  that  tables  of  Z-transforms  may  be  de- 
veloped in  the  same  manner  as  tables  of  ordinary  Laplace  transforms. 

One  of  the  big  advantages  of  the  Z-transform  is  that  it  is  possible 
to  obtain  the  pulsed  output  of  a  linear  system  by  expressing  a  pulse  trans- 
fer function  which  relates  the  input  and  output  Z-transforms.   The  pulse 
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transfer  functions  can  be  expressed  in  closed  form  and  tables  of  Z-trans- 
forms  are  available  in  most  standard  texts  covering  sampled-data  systems. 
An  abbreviated  table  is  given  in  Table  6-1. 
6.3  Application  to  Linear  Continuous  Systems 

The  application  of  the  numerical  method  using  Z- transform 
technique  will  be  illustrated  in  a  linear  continuous  feedback  system  as 
shown  in  figure  6-2(a).   The  system  is  mathematically  converted  into  a 
sampled  system  by  placing  synchronized  samplers  at  input  and  output  and 
at  point  X  in  the  feedback  path  in  figure  6-2 (b).   The  sampling  frequency 
should  be  high  enough  to  cause  the  error  introduced  by  the  sampled  system 
to  be  small.   By  locating  the  sampler  in  the  feedback  loop  and  requiring 
that  the  transfer  function  G(s)  be  low  pass  in  nature,  then  only  the  for- 
ward gain  of  the  original  system  need  be  considered.   Later  examples  will 
show  the  effect  of  the  ratio  of  sampling  frequency  uJ   to  the  Bode  cross- 
over frequency  ul>c  of  the  original  system. 

Following  the  sampler  S  in  the  feedback  loop  is  the  element  B(s)j 
whose  function  is  to  reconstruct  the  original  output  of  the  forward  path 
as  faithfully  as  possible.   This  element  reconstructs  a  polygonal  approxi- 
mation to  the  actual  function  by  converting  each  pulse  into  a  triangular 
generating  function  as  shown  in  figure  6-3.   This  requires  the  response  of 
B  to  a  short  pulse  or  impulse  of  the  same  area  to  be  a  triangle  as  shown 
in  figure  6-4(a).   This  response  is  not  physically  attainable  since  it 
requires  a  slope  at  negative  time,  but  this  has  no  effect  on  the  purely 
mathematical  solution. 

The  combination  of  sampling  switch  operating  with  a  period  T 
and  the  triangle  function  generator  B  replace  the  direct  connection 
which  existed  in  the  original  system.  The  errors  introduced  into  the 
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.    Table  6-/ 
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Figure  6-2.     (a)     Continuous  Feedback  System 


Sampled  Model  for  Computation 
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•(t) 


Polygonal 
Approximation 


Generating 
Triangle 


Figure  6-3.     Polygonal  Approximation  of  a  Function  f (t) 
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Figure  6-4..   (a)  Generating  Triangle 

(b)  Ramps  Which  Produce  Triangle 
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solution  are  a  result  of  the  Inability  of  this  combination  to  faithfully 
reproduce  the  feedback  signal. 

The  transfer  function  B(s)  of  the  triangle  function  generator  is 
obtained  from  the  Laplace  transform  of  the  triangle  function  in  figure 
6-4(a).   This  triangle  is  a  result  of  the  combination  of  three  ramp  func- 
tions as  shown  in  figure  6-4(b).   The  Laplace  transform  of  these  three 
functions  is:  ,    t-c     ^  \  -Ts 


-  -  _g(l-e-r 


B(s^  &X(\-z 


6.3 


The  relation  between  the  Z-transform  of  the  output  and  input  in 
figure  6-2 (b)  is  given  by 

©*(Z)   =      &&*(*) 6.4 

and  yields  the  pulse  transfer  functions  for  which  tables  are  available. 

If  the  response  of  a  linear  system  to  an  input  F(s)  is  desired, 
first  replace  it  with  the  system  shown  in  figure  6-2 (b).   A  proper  choice 
of  sampling  frequency  (x^  is  made.   The  pulse  transfer  functions  RG*(z) 
and  HBG*(z)  are  computed,  and  combined  according  to  the  basic  laws  of 
algebra,  yielding  a  form  of  a  ratio  of  polynomials  in  powers  of  z.   To 
obtain  the  output  pulse  sequence  in  the  time  domain,  a  simple  process  of 
long  division  is  required.   The  division  process  will  result  in  a  power 
series  of  the  form 

A,(2)=  <x0*  alzi+aL£z+  a3^+aM2Ti  ■ .  »  +  an2'V—    6.5 

The  magnitude  of  the  output  in  the  time  domain  is  just  the  coefficient 

of  each  of  the  terms  in  this  series  at  a  time  corresponding  to  the 

power  of  z.   Thus  the  output  at  the  fourth  sampling  instant  is  a,  in  the 

expansion. 
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The  long  division  process  is  continued  until  the  desired 

number  of  points  have  been  obtained.   The  effect  of  the  sampling  period 

T  ,  or  of  the  sampling  frequency  uJ<   will  be  elaborated  on  in  section  6.5. 
s 

6.4  Application  to  Linear  Continuous  Systems  with  Transportation  Lag. 

It  has  been  shown  in  section  4.1  that  the  transfer  function  of  a 

-Ts 
transportation  lag  or  simple  dead  time  lag  is  merely  e     where  T  is  the 

value  Tp  ,  or  the  time  lag. 

If  this  transfer  function  is  inserted  into  the  linear  system  of 

figure  6-2(b)  as  G. (s)  then  the  solution  for  the  output  will  become 

The  value  of  this  numerical  method  employing  pulse  transforms 

-T  s 
is  now  easily  seen.  The  z-transform  of  the  time  delay  e   °   simply  be- 
comes z   from  the  basic  definition  of  the  z-transform  itself.   Since  the 
remaining  terms  of    t:*(z)  are  also  polynomials  in  powers  of  z,  this 
does  not  add  any  more  complexity  to  the  basic  long  division  solution  of  the 
transient  response. 

It  must  be  noted  that  the  T  introduced  in  the  polynomial  recon- 
struction and  the  z-transforms  in  Table  6-1  is  T  ,  the  sampling  interval. 

9 

This  must  be  distinguished  and  kept  separate  from  T  ^  ,the  time  delay. 

One  of  the  basic  rules  for  the  method  is  that  the  sampling  time 
T  must  be  equal  to  or  an  integral  sub-multiple  of  the  time  delay  T  ~  . 
It  will  be  necessary  to  express  6)*(z)  only  in  terms  of  the  sampling  time 
T  in  order  that  the  value  of  the  transient  response  will  be  obtained  at 
each  sampling  interval.  This  method  will  be  shown  in  the  sample  solu- 
tions in  section  7. 
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6.5   Requirements  for  Reasonable  Accuracy 

The  ratio  of  U^  to  the  crossover  frequency  uJQ    of  the  open 

loop  transfer  function  will  determine  the  relative  accuracy  of  the  system. 

The  comparison  of  sampling  period  T  with  the  expected  settling  time  and/ 

s 

or  time  delay  in  the  system  will  determine  whether  the  method  will  supply 
the  required  answer  within  a  reasonable  number  of  long  division  steps. 

Ideally,  to  obtain  extremely  good  accuracy  (within  17.  of  known 
answer)  it  is  necessary  to  select  a  value  of  u^  which  is  twice  the  fre- 
quency at  which  the  open  loop  system  is  down  40  db  on  the  Bode  diagram. 
For  the  systems  investigated  here  and  explained  in  section  2,  u>0  is  approxi- 
mately 1  rad/sec.   For  the  system  where  G(s)  =  —  ,  this  would  require  u-^ 

s 

=  20*0  rad/sec.   Since  T  =  —   this  would  require  a  T  of  approximately 

0.03  sec.   With  a  normal  settling  time  for  this  system  of  4  seconds  this 

4 
would  require  a  total  of  —   or  133  long  division  steps  to  obtain  the 

1  UJ 

transient  value  at  the  end  of  4  seconds.   In  order  to  reduce  the  computa- 
tional work  and  to  remain  within  approximately  an  accuracy  of  5%,  vari- 
ous ratios  of  u)     to  toc    were  selected  and  results  tabulated  against  re- 
sults of  former  traditional  methods  of  solution. 

The  problem  of  reducing  computational  work  is  further  aggravated 
by  the  introduction  of  the  time  delay  into  the  system.   One  basic  re- 
quirement of  this  numerical  method  is  that  the  sampling  time  T   must  be 

s 

equal  to  or  an  integral  sub-multiple  of  the  time  delay  T  0      .   For  systems 

with  a  long  settling  time,  for  example  5  seconds,  and  a  short  time  delay 

of  0.10  seconds,  it  would  require  50  long  division  steps  to  produce  the 

value  of  the  transient  response  at  5  seconds.   This  is  only  if  it  would 

be  possible  to  select  T  =  T  -   .   In  the  event  that  the  value  of  T  must 
r  s     0  s 

be  chosen  smaller  to  satisfy  the  requirement  of  the  ratio   u^   to  u>0   for 
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accuracy,  then  it  would  require  50n  long  division  steps,  where 

To 


n  = 


T, 


It  will  be  demonstrated  in  the  next  section  that  this  method 
will  provide  reasonable  engineering  accuracy  (approximately  57.)  with  a 
small  investment  of  time.   This  can  be  accomplished  with  a  ratio  of  tO^ 

to  n>0    as  low  as  possible,  and  a  T  which  is  the  lowest  sub-multiple  of 

™        ■>  2n 

T     as  long  as  xt    is  not  greater  than  r—  . 

This  method  can  give  extremely  good  accuracy  at  the  expense  of 
an  extremely  long  and  repetitious  series  of  calculations.   Its  advantage 
lies  in  the  fact  that  where  the  ratios  mentioned  above  are  favorable,  it 
should  be  possible  to  compute  the  transient  response  of  a  system  in  a 
short  period  of  time  with  no  more  powerful  tools  than  algebra  and  a  slide 
rule  or  desk  calculator. 

The  next  section  will  demonstrate  the  accuracy  and  simplicity 
of  the  method  in  solving  simple  feedback  systems.   Computational  time 
and  labor  will  be  compared  with  other  methods  currently  available  to  the 
engineer. 
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7.    Comparison  of  Numerical  Method  with  Other  Methods  of  Solution 
7. 1   Numerical  Method  Solution 

The  numerical  method  solution  to  the  second  order  system  des- 
cribed in  section  2  will  now  be  demonstrated  in  detail.   The  system  to 
be  solved  is  shown  in  figure  6-2 (b).   The  input  is  a  unit  step.   The  time 
delay  will  be  one  second.   The  sampling  time  interval  will  also  be  one 
second.   The  table  of  results  in  section  7.5  will  also  show  the  solution 
results  for  other  sampling  intervals. 

The  output  Qq.O)      in  the  Laplace  notation  is 

to  e'?TD 


^  ^  J  /ocl-*t>BCs) 


7.1 


/* 


5(S+'o) 


where 


e-5    /  -st|2 


Bu)  =     -?-(,-&") 


Ts>  v        ^      '  7-2 


Combining  -sX, 

toe 


r\    .  .     _       Sz  (s-h/O) 

£<:(*)     ~  .  /q(/-£-^  7.3 

'    r        S3  (S+/o) 
It    is   important  to  note  that   the  T  in  the  time  delay  represents 

TD     ,    the  delay  time,    whereas   the  T   in  the  expression   for  B(s)   and  also   in 

the   Z  transforms   represents  T   ,    the   sampling  time.      In  the  algebraic  work 

s 

in  combining  the  expressions  it  will  be  necessary  to  express  T&   in  terms 

of  T   and  obtain  the  output  equation  in  terms  of  T   only.   In  this  ex- 

s  s 

ample,  since  To   equals  T   equals  one  second,  it  will  not  be  necessary  to 

s 

account  for  this  difference. 

By  expressing  numerator  and  denominator  separately  in  partial 
fraction  expansions,  and  then  transforming  each  separately  by  means  of 
the  transforms  in  table  I  of  section  6  we  obtain 
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ec(2)  =  ioo-^-r 7.4 

zoo(i~z-') 
Combining  numerator  and  denominator  results    in 


0    (£)    -        iW'+ZOr3 7.5 


By  carrying  out  the  long  division  indicated  in  equation  7.5 
we  obtain  a  power  series  in  Z 

6c  (SO  =  0.12*+ M3t 2*+  i.385-£~'¥  +  0.9/  Z?5^  •  •  •        7>6 

We  may  directly  transform  this  expression  into  the  time  domain. 
It  is  not  really  necessary  since  the  only  effect  on  the  output  will  be 
to  express  the  series  in  terms  of  the  sampling  periods.   In  either  case 
the  coefficients  are  the  same  and  are  the  value  of  the  transient  response 
at  each  sampling  instant. 

Qc(t)   -  M(fc-2T)+/.53/(t-JT-)+  J.38S(t-HT)+   •  •  •     7.7 

The  introduction  of  a  different  time  delay  or  sampling  time 

into  the  problem  merely  changes  the  algebraic  relationship  between  TD 

and  T   and  once  this  accounting  has  been  accomplished,  the  algebra  and 
s 

final  solution  work  is  similar  to  that  shown  here.   No  additional  diffi- 
culty is  introduced  into  the  problem. 
7.2  Power  Series 
7.21  First  Order  System 

The  solution  of  the  transient  response  of  a  first  order  system 
with  transportation  lag  by  the  employment  of  a  power  series  expansion  is 
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quite  straight  forward,  although  tedious  and  time  consuming.   The 
system  under  consideration  is  that  shown  in  figure  7- 1(a)  where  T  is 
equal  to  one  second. 

The  closed  loop  transfer  function  may  be  written 


&c 

= 

S 

Q* 

/  +  £— 

For  a 

unit 

step  input 

ft 

/ 

7.8 


7.9 


Then 


~T5    /  ^-T5\-/ 


ft   ,    -3L-       =      e     //  +  _e_J) 


7.10 


f+S 


-Ts    -1 
The  quantity  (  1  +  e    )       may  be  expanded  in  a  power  series 


with  the  following  result. 

-Ts 


ec  = 


_   e 


2. 


G-Ts        c-iTs       C-3TS       e-vr,  -j 

/-  —  +-^r  --55-  +  -J3  j        7.11 


n        e~r5     C'zTs     e'3Ts     c-vTi      e"575 

be  -    -jt        —^3-  +  — ^ ^F"  -* y$-  —  ■  -  •        7oi2 


Since  T  =   1 


tfc-      ^  S3   r      Sh  ^  s6 

Taking  the  inverse  Laplace  transform  of  each  term  yields 
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7.13 


7.14 


Figure  7-l(a) .  First  Order  Feedback  Control  System 
V.'ith  Transportation  Lag. 


Figure  7-1  (b). 


Second  Order  Feedback  Control  System 
With  Transportation  Lag. 
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The  transient  response  may  now  be  determined  by  substituting 
values  of  t  into  equation  7.14   .   The  magnitude  of  the  output  at  any 
time  t  may  be  determined  as  exactly  as  desired  by  evaluating  all  terms 
and  by  evaluating  each  term  to  a  sufficient  number  of  decimal  places. 
The  amount  of  time  required  to  obtain  a  transient  solution  by  this  method 
increases  directly  as  the  time  interval  between  sampling  instants  is  de- 
creased.  It  also  increases  as  the  settling  time  increases  since  increased 
settling  time  increases  the  number  of  terms  which  must  be  evaluated  at 
higher  values  of  t.   The  sample  problems  solved  by  this  method  required 
about  twelve  hours  to  complete. 
7.22  Second  Order  System 

While  the  use  of  a  power  series  expansion  is  a  straight  forward 
method  of  determining  the  transient  response  of  a  first  order  system  with 
transportation  lag,  the  method  requires  far  too  much  work  to  be  considered 
practical  for  higher  order  systems.   The  following  example  will  demons- 
trate why  the  method  fails  as  a  practical  working  tool.   The  system  under 
consideration  is  that  shown  in  figure  7-l(b)  where  T  is  equal  to  one 
second. 


The  closed  loop  transfer  function  may  be  written 
I  + 


Or  '  i  _.       er5T"  7a5 


For  a  unit   step   input 


s 


Then  cT 

e 


sT      /  _ST    \~i 


Uc  "    HFF—  -       ^/ ,      /+  TZTT-Tx  7.16 


'         S(o. 


iS+\) 


46 


The  quantity   /+  — —       may  be  expanded  in  a 

power  series  with  the  following  results. 


c-isT  ^-H5T  7.17 

V(a/s+o3  +    ~E*(ais±i) 

At  this  point  it  may  readily  be  seen  that  the  task  of  determin- 
ing the  inverse  Laplace  transform  of  each  term  is  momentus  and  the  work 
involved  increases  tremendously  with  each  increase  in  the  order  of  the 
term.   Thus  it  is  maintained  that  the  use  of  a  power  series  is  impracti- 
cal in  determining  the  transient  response  of  systems  of  second  order  or 
higher  which  are  complicated  by  transportation  lag. 
7. 3  Pade  Approximation  (Arithmetic) 

The  transient  response  will  now  be  determined  for  the  system 
of  section  7.1  by  simulating  the  Laplace  shift  operator  with  the  second 
order  Pade  approximant.   The  major  steps  in  the  solution  will  be  outlined 
but  the  tedious  arithmetical  calculations  will  of  necessity  not  be  repeated 
in  their  entirety. 

The  system  in  block  diagram  form  is  shown  in  Figure  7- 1(b). 

Since  the  duration  of  the  transportation  lag  is  one  second  the  Laplace 

—  s 
shift  operator  becomes  e   .   As  indicated  in  section  4.22 

e"5^      s2-^  +  >2~ 7.i8 

Substituting  the  Pade  approximant  in  for  the  Laplace  shift 
operator  and  manipulating  by  block  diagram  algebra  yields  the  closed 
loop  transfer  function  for  the  system, 

-^-     -  S(S-Ho)(5N-6,S-H2. 

9R  /    +  /0(S*-4s  +  /2) 

/  S(5-HoXsV6s+iO 
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For  a  step  Input,  Q    -     — 


n    -  J.    /    lOS^-LOS  +  tZQ )  .  7n 

The  transient  response  of  this  system  to  a  step  input  is  con- 
tained in  Laplace  notation  above.  The  problem  of  reverting  to  the  time 
domain  remains. 

The  first  step  is  the  determination  of  the  four  roots  of 

4      3      2 
S  +  16S  +  82S  +  60S  +  120.   Since  these  roots  are  two  complex  conju- 
gate pairs  whose  real  and  imaginary  parts  are  not  necessarily  integral 

the  Lin  method  was  selected  for  their  solution.   The  first  trial  divisor  is 

2  2 

82S  +  60S  +  120  which  equals  S  +  .733S  +  1.46.   Five  iterations  were 

2  2 

necessary  to  yield  factors  of  (  S  +  .463S  +  1.65  )  (  S  +  15.548^+  73.2  ). 

Application  of  the  quadratic  formula  further  reduced  the  factors  to  (  S  + 
.232  +  jl.264  )  (  S  +  7.77  +  J3.59  ). 

Next,  the  method  of  residues  was  employed  to  determine  the  co- 
efficient of  the  inverse  transform  of  each  of  the  five  poles. 
Thus 


9C  (t)  =  /  ♦  (-j.*  H.m)c(-^ii!cv)c,  (,„,^,)  er-«>  W 


7.21 


This   further  reduces  to 


-JSIC     7?1Z(Z  cos3sn  +  2.6J  sin3.SU)  7'22 


The  full  transient  response  of  the  system  may  now  be  plotted  by 
evaluating  Qc(t-)   at  various  values  of  t. 
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The  time  required  for  the  solution  of  the  transient  response 
of  the  system  by  this  method  is  between  seven  and  eight  hours.   It  is 
worthy  of  note  that  the  roots  of  the  characteristic  equation  determined 
by  this  method  are  reasonably  close  to  the  values  obtained  from  the  more 
accurate  but  more  lengthy  root  locus  solution.   Arithmetic  solutions 
making  use  of  the  Fade  approximation  were  limited  to  the  use  of  the 
second  and  third  order  approximants  because  the  time  required  to  solve 
a  problem  when  higher  order  approximants  are  used  becomes  unreasonably 
long. 
7.4  Pade  Approximation  (Computer) 

In  section  4.22  it  was  pointed  out  that  the  fourth  order  Pade 
approximant  was  more  than  adequate  for  handling  the  vast  majority  of  prob- 
lems involving  transportation  lag.   Accordingly,  all  problems  solved  by 

arithmetical  techniques  were  also  solved  on  the  Donner  Analog  Computer 

-sT 
using  the  fourth  order  Pade  approximant  to  simulate  e  '  .   The  solution 

to  the  problem  by  this  method  was  established  as  the  control  against 
which  the  arithmetic  solutions  were  compared.   The  fourth  order  approxi- 
mant was  chosen  because  it  occupied  only  six  amplifiers  on  the  Donner 
Computer  thus  leaving  four  amplifiers  free  on  which  the  remainder  of  the 
system  could  be  programmed.   The  computer  program  for  the  fourth  order 
Pade  approximant  is  contained  in  the  Appendix.   This  program  was  specifi- 
cally written  for  a  range  of  values  of  transportation  lag  between  0.1 
seconds  to  1.0  seconds.   In  each  amplifier  only  one  resistor  R T   had  to 
be  changed  to  change  the  value  of  the  transportation  lag  between  the  limits 
stated  above.   For  transportation  lags  outside  of  these  limits  the  program 
should  be  modified  so  that  EL,  does  not  exceed  the  bounds  of  0.1  to  1.0 

megohm. 
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7.5  Table  of  Results 

This  table  of  results  is  an  intent  to  show  the  relative  ac- 
curacy of  the  various  methods,  and  to  compare  the  solution  time  in- 
volved to  achieve  this  accuracy.   The  first  order  system  described  in 
section  2  was  used  to  establish  the  basic  comparisons  for  accuracy,  and 
to  demonstrate  the  feasibility  of  the  numerical  method.   Results  indi- 
cate that  the  numerical  method  approach  provides  results  of  comparable 
accuracy  to  other  methods  with  substantial  savings  in  computational  time. 
Results  also  indicate  that  increased  accuracy  may  be  obtained  with  the 
numerical  methods  solution  by  increasing  the  sampling  frequency.   It 
must  be  pointed  out,  however,  that  the  time  to  arrive  at  the  solution 
increases  directly  in  proportion  to  the  increase  in  the  sampling  fre- 
quency. 

Table  7-1  contains  the  basic  information  from  which  general 
conclusions  about  relative  accuracy  will  be  based.   The  results  indicate 
the  value  of  the  transient  response  for  a  first  order  system  with  a  one 
second  time  delay  as  obtained  by  various  methods.   Results  are  shown  at 
one  second  intervals.   The  power  series  solution  is  the  most  accurate 
available  and  will  be  taken  as  the  standard  by  which  the  other  methods 
will  be  compared. 

Table  7-2  contains  similar  information  for  the  first  order  sys- 
tem with  a  0.5  second  delay.   The  numerical  method  solution  for  a  ratio 
of  (jjs  to  <-0o  ratio  of  12.56  is  compared  because  it  achieves  the  stated 
requirements  of  the  study  of  less  than  57.  deviation.   It  will  be  the  ratio 
used  in  higher  order  systems. 

Tables  7-3  and  7-4  display  results  for  the  second  and  third  order 
systems  respectively  for  a  0.5  second  delay.   The  basis  for  comparison  in 
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these  systems  was  the  computer  solution  employing  the  fourth  order  Pade 
approximant.   The  deviations  here  from  the  standard  are  of  the  same 
orders  of  magnitude  as  experienced  in  the  first  order  system.   In  these 
results,  a  plus  deviation  indicates  a  result  higher  than  the  standard, 
and  a  minus  deviation  a  result  lower.   The  numerical  method  solution 
when  compared  with  the  standard  characteristically  maintains  the  same 
sign  at  the  same  relative  point  of  the  transient  response.   That  is,  the 
transient  response  as  determined  by  the  numerical  method  coincides  with 
or  leads  the  standard,  and  would  therefore  tend  to  show  a  slightly  short- 
er rise  time  and  time  to  peak  overshoot.   Its  accuracy  in  determining  the 
peak  overshoot  and  settling  time  was  well  within  the  engineering  accuracy 
prescribed  for  this  investigation. 

It  must  be  understood  that  the  solutions  by  the  numerical  methods 
technique  yield  "exact"  values  only  at  the  sampling  instants.   Interpola- 
tion may  be  done  to  determine  information  at  times  between  these  sampling 
instants  with  assurance  of  reasonable  accuracy  as  the  magnitude  and  alter- 
nating signs  of  the  deviations  of  the  tabulated  results  from  the  standard 
indicate. 

All  solutions,  both  computer  and  analytic,  employing  the  Pade 
approximation  exhibit  large  deviations  at  the  first  few  tabulated  times,. 
These  are  due  to  the  "ringing"  discussed  in  section  4.22  and  illustrated 

in  figure  A- 3.   This  is  an  unavoidable  disadvantage  of  the  Pade  approxi- 

-  -sT 
mat ion  of  e 

It  cannot  be  over  emphasized  that  the  time  required  to  obtain 

the  transient  solution  by  the  numerical  method  was  in  most  cases  from 

one-fifth  to  one-tenth  the  time  required  by  previously  used  analytical 
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and  graphical  methods.   Its  value  also  lies  in  the  fact  that  the  pro- 
cedure is  simple  and  easy  to  follow.   Once  the  engineer  or  student 
becomes  familiar  with  the  procedures,  it  is  possible  to  recognize  and 
correct  computational  errors  before  they  are  compounded  in  following 
steps.   This  is  simply  not  possible  in  the  other  analytical  methods 
until  the  final  step  of  the  solution.    Similarly,  a  change  in  the 
transportation  lag  or  the  sampling  interval  is  easily  handled,  and 
a  great  deal  of  the  algebraic  manipulations  need  be  done  only  once 
for  each  basic  system. 
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TABLE   7-1 
COMPARISON  OF  TRANSIENT    RESPONSE  AS   OBTAINED  BY 
POWER  SERIES   WITH  OTHER  ARITHMETIC  METHODS 
FIRST  ORDER   SYSTEM  -    TRANSPORTATION   LAG  EQUAL  TO  ONE   SECOND 


Time     Power 
Sec.    Series 


2nd 
Order 

fArlth.) 


7«  Num. 

Dev.      Meth. 
ojs1  6.28 


7.                Num.               7o            Num.               7. 
Dev          ,    Meth.          Dev.    .     Meth.        Dev. 
u-fe-12.56     oj0-25.12      


0 

0 

.011 

- 

0 

- 

0 

- 

0 

- 

1 

0 

.291 

- 

0 

- 

0 

- 

0 

- 

2 

1.00 

1.013 

+1.3 

1.00 

0 

1.00 

0 

1.00 

0 

3 

1.50 

1.485 

-1.0 

1.50 

0 

1.50 

0 

1.50 

0 

4 

1.167 

1.158 

-0.8 

1.250 

+7.1 

1.1875 

+1.8 

1.1718 

+0.41 

5 

.791 

.795 

+0.5 

.875 

+10.6 

.8125 

+2.7 

.7969 

+0.75 

6 

.8416 

.848 

+0.8 

.8125 

-3.5 

.8320 

-1.1 

.8392 

-0.28 

7 

1.057 

1.055 

-0.2 

.9688 

-8.3 

1.0331 

-2.3 

1.0509 

-0.58 

8 

1.100 

1.098 

-0.2 

1.0782 

-2.0 

1.0993 

-0.1 

1.1022 

+0.2 

9 

1.004 

1.005 

+0.1 

1.0548 

+4.9 

1.0205 

+1.6 

1.0087 

+0.37 

10 

.948 

.950 

+0.2 

.9884 

+4.2 

.9564 

+0.9 

.9543 

+0.66 

11 

.9798 

.981 

+0.1 

.9669 

-1.3 

.9719 

-0.8 

.9916 

+1,20 

12 

1.021 

1.019 

-0.2 

.9894 

-3.1 

1.0108 

-1.0 

1.0329 

+  1.15 
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TABLE   7-2 
COMPARISON  OF  TRANSIENT   RESPONSE  AS   OBTAINED 
BY  POWER  SERIES   WITH  NUMERICAL  METHOD   RESULT 
FIRST  ORDER  SYSTEM-TRANSPORTATION   LAG  EQUAL  TO  0.5   SECOND 


Time  Power  Numerical                                   % 

Method 

Sec  Series  Dev. 

^ -12.56 

0  0  0                                        0 

0.5  0  0                                        0 

1.0  0.5  0.5                                    0 

1.5  .875  .875                                 0 

2.0  1.028  1.031  +0.29 

2.5  1.039  1.055  +1.06 

3.0  1.021  1.033  +1.18 

3.5  1.006  1.011  +0.5 

4.0  .999  1.000  +0.1 
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TABLE   7-3 
COMPARISON  OF  TRANSIENT    RESPONSE   FROM  COMPUTER  SOLUTION 
WITH  NUMERICAL  METHOD  AND  PADE  APPROXIMANT    (ARITHMETIC) 
SECOND  ORDER  SYSTEM  -    TRANSPORTATION   LAG  EQUAL  TO  0.5    SECOND 


Time 
Sec 

Computer 

(jOo 

Numerical 
Method 

^12.56 

% 
Dev 

Pade 
Approx. 
(Arith) 

% 
Dev. 

0 

0 

0 

- 

.01 

- 

0.5 

0 

0 

- 

.0018 

- 

1.0 

.31 

.35 

+12.8 

.404 

+30.4 

1.5 

.77 

.803 

4  4.3 

.846 

+  9.9 

2.0 

1.04 

1.078 

+  3.7 

1.073 

+  3.2 

2.5 

1.115 

1.153 

+  3.4 

1.111 

-  0.5 

3.0 

1.100 

1.110 

+  0.9 

1.066 

-  3.1 

3.5 

1.050 

1.037 

-  1.2 

1.017 

-  3.1 

4.0 

1.015 

.984 

-  3.1 

.992 

-  2.3 

55 


TABLE   7-4 
COMPARISON  OF  TRANSIENT   RESPONSE    FROM  COMPUTER  SOLUTION 
WITH  NUMERICAL  METHOD  AND  PADE  APPROXIMANT    (ARITHMETIC) 
THIRD  ORDER  SYSTEM   -    TRANSPORTATION   LAG  EQUAL  TO  0.5   SECOND 


Time 

Numerical 

Pade 
Approx. 
(Arith) 

% 
Dev. 

Sec. 

Computer 

Method 
■=12.56 

7. 
Dev. 

0 

0 

0 

- 

.146 

- 

0.5 

0 

0 

- 

.0203 

- 

1.0 

.26 

.35 

+34.6 

.359 

+38.0 

1.5 

.73 

.803 

+  10.0 

.821 

+12.5 

2.0 

1.050 

1.078 

+  2.7 

1.09 

+  3.8 

2.5 

1.150 

1.153 

+  0.3 

1.150 

0 

3.0 

1.115 

1.110 

-  0.5 

1.099 

-  1.4 

3.5 

1.060 

1.037 

-  2.2 

1.029 

-  2.9 

4.0 

1.005 

.984 

-  2.1 

.987 

-  1.8 
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8.    Conclusions  and  Recommendations  for  Further  Investigation 

The  investigation  into  the  effect  of  transportation  lag  on 
feedback  control  systems  reported  on  in  this  thesis  has  demonstrated 
four  ideas  of  practical  significance.   It  has  also  opened  up  a  number 
of  avenues  for  further  investigation. 

The  primary  conclusion  which  the  investigation  supported  is 
the  applicability  of  the  numerical  method  technique  devised  by  RAGAZZINI 
and  BERGEN  to  the  problem  of  determining  the  transient  response  of  linear 
feedback  control  systems  which  are  complicated  by  transportation  lag.   Of 
all  the  mathematical  and  graphical  tools  available  to  the  engineer  it 
offers  the  optimum  combination  of  speed,  accuracy,  and  simplicity.   It 
may  be  used  by  students  whose  background  includes  courses  only  as  advanc- 
ed as  linear  feedback  control  systems,  provided  the  Z-transforra  is  included, 

The  deviation  of  the  transient  response  as  determined  by  the 
numerical  method  from  that  determined  by  the  techniques  currently  in  use 
may  be  held  to  less  than  three  percent.   Meanwhile  the  time  required  to 
determine  the  transient  response  has  been  reduced  by  more  than  seventy- 
five  percent.   Thus  the  decision  of  whether  compensation  is  required  or 
not  may  be  made  much  more  rapidly  than  before.   The  form  the  compensation 
will  take  will,  of  course,  vary  from  problem  to  problem  and  from  designer 
to  designer. 

The  suitability  of  the  numerical  method  to  the  particular  pro- 
blem of  transportation  lag  suggests  an  area  of  further  investigation.   In 
determining  the  transient  response  the  system  has  been  mathematically  con- 
verted to  a  sampled-data  system,  and  is  described  by  Z-transforms.   In 
order  to  determine  whether  the  system  will  meet  specified  performance 
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criteria  it  is  necessary  to  calculate  the  transient  response.   If  the 
system  is  not  satisfactory,   compensation  techniques  applicable  to  con- 
tinuous feedback  systems  would  be  employed  and  the  transient  response  of 
the  compensated  system  recalculated.   However,  since  the  system  is  al- 
ready a  mathematical  sampled-data  system,  it  would  appear  that  the  suit- 
ability of  sampled-data  system  performance  criteria  and  compensation 
techniques  to  the  problem  is  worthy  of  investigation.   Should  these 
techniques  prove  successful  the  reverting  back  and  forth  from  the  Laplace 
transform  to  the  Z- trans  form  for  compensation  could  be  avoided  and  a 
greater  saving  in  time  and  tedium  effected. 

The  second  noteworthy  conclusion  of  this  investigation  is 
that  the  Bode  diagram  not  only  offers  the  fastest  tool  for  the  deter- 
mination of  system  stability  but  it  also  will  indicate  which  set  of  roots 
will  go  unstable  first.  This  information  is  made  available  by  no  other 
technique  except  the  root  locus  and  a  conservative  estimate  would  place 
the  time  differential  using  the  root  locus  at  a  full  order  of  magnitude 
higher  than  the  Bode. 

The  criterion  of  dominance  of  a  set  of  roots  in  a  higher  order 
system  provides  the  basis  for  the  extension  of  the  performance  character- 
istics of  a  second  order  system  to  these  higher  order  systems.   From  the 
information  in  the  Bode  and  root  locus  graphical  plots,  it  may  be  practical 
to  construct  a  set  of  second  order  curves  with  transportation  lag  which 
might  have  applicability  to  the  analysis  of  higher  order  systems, 

A  third  significant  conclusion  is  that  the  Pade  approximation 
for  the  Laplace  shift  operator  is  considerably  more  linear  in  its  phase 
shift  characteristics  than  the  recent  analog  computer  literature  would 
indicate.   An  inspection  of  figure  4-4  will  reveal  that  the  Pade  approxi- 
mation provides  a  completely  linear  phase  shift  until  that  point  where  it 

58 


deviates  from  the  ideal.   Thus,  the  simulation  of  time  delay  on  the  analog 
computer  requires  merely  the  choice  of  a  high  enough  order  Pade  approxi- 
mant  which  will  linearly  phase  shift  the  frequencies  desired.   It  would 
appear  that  the  curve  fitting  techniques  currently  in  existence  are  not 
justified  due  to  their  increased  complexity  and  lack  of  a  simple  repeti- 
tive scheme  for  programming, 

A  fourth  interesting  result  of  this  study  is  the  information  on 
stability  limits  for  the  second  order  system.   When  plotted  on  logarithmic 
cross-section  paper  as  shown  in  Figure  5-4,  the  non-dimensionalized  gain 
vs.  the  time  delay  produced  a  straight  line  curve  from  which  it  was  possible 
to  write  an  equation  describing  the  relationship  between  gain  at  stability 
limit  and  transportation  lag  for  any  second  order  system. 

This  result  is  obtainable  from  root  locus  or  frequency  response  plots 
and  is  easily  presented  because  of  the  ease  in  non-dimensionalizing  a 
second  order  system.   Here  again,  a  further  extension  of  this  informa- 
tion would  be  possible  for  higher  order  systems.   The  method  of  display- 
ing the  information  would  depend  on  the  way  in  which  the  additional  time 
constants  of  the  higher  order  systems  were  handled. 
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9.   Appendix  -  Computer  Program  for  Fourth  Order  Pade  Approximant 

The  Pade  approximation  of  e  is  discussed  in  Chapter  4.   An 

analog  computer  program  for  the  fourth  order  approximant  which  was  used 

in  the  study  is  developed  here. 

The  Pade  approximation  states  that 

C    T   /u  +  \^-»oo    ^      ,  .  ?'1 


where 


(**„)—-  &  x) 


r   /v\  .  1  +    vx        v(v-i)  xz 
r^     '      (uw^ll      (u*v)(u*v-.)2! 


.  ..  -r 


v(v-i)  —  a-i  x* 

(U+V)(U+V-|)  •  ••  (lA-H)Vj 


and 


.  zo        izo  x^o  i(*%o 

£-  sT  -     /  -  TT  "*"  7sT)2  ~  (ST)3  "*"   (ST}' 


sT  +  STY-  +  (st)3  +  Tiry 


As  shown  in  Section  4.22 


e-sT=     jT(f(t-T)] 


9.2 


u'v  '      (v+n.)f.'      (v+u)(m-i)ai 

...  +(_j)*  u(u-0-— e-ix* 

Substituting  u  =  v  =  4,  and  x  =  -sT  yields 

-sT=    (ST)4-  Z0(ST)*+  /gOftT)*-  g*/0(sT)  +  /^0       ^ 
(ST}Vao(sT)3+/?0(sT)24^0(5T)+/&&0 


9.5 


9.6 
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Since  the  differential  operator  and  Laplace  operator  are  inter- 
changeable when  initial  conditions  are  zero,  equations  9.5  and  9.6  may  be 
combined  as  follows 


f(t-T) 


/-  -2o 


ISO 


Sf-O 


-£ill__£LH 


/4>8Q 
DHT+ 


I 


20 


/go 


D3T3 


9.7 


The  final  form  of  equation  9.7  from  which  the  analog  computer  program 
shown  in  figure  9-1  was  developed  is 


+ 


f  Kt)-iff(t-r)]i 


-[^4(t)4^(t-T)]^ 


D 


k?o 


l(*20 


-rt  1  C^-3   -j-  4- 


{(t-T) 


D' 


9.8 
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TABLE  OF  COMPUTER  COMPONENTS 
Potentiometer  Settings  Component  Values 

a.  »  Rici  1  Megohm 

a3  "  R3C2 
a4  =  0.7  C2 

a5  -  R5C2 

a6  -  R6C3 
a  -  0.467  C 

38  "  R8C3 

a9  B  R9C4 

310  "  °-9  °4  Rf5 

all  *  R11C4 

1  9 

12   - —  1  Microfarad 

Rf5 

a13=i^  Cl    C2 

al4  * ■  !tt  C3    C4 

Rf6 


The  time  is  scaled  such  that  computer  time  is  equal 
to  ten  times  problem  time. 


Rl 

R2 

R3 

R5 

R6 

R8 

R9 

Rll 

Ru 

R14 

Rf6 

2  Me 

gohm 
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